Centrality bin size dependence of multiplicity correlation 
in central Au+Au collisions at ^snn=200 GeV 
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We have studied the centrahty bin size dependence of charged particle forward-backward multi- 
plicity correlation strength in 5%, 0-5%, and 0-10% most central Au-|-Au collisions at ysNN=200 
GeV with a parton and hadron cascade model, PACIAE based on PYTHIA. The real (total), sta- 
tistical, and NBD (Negative Binomial Distribution) correlation strengths are calculated by the real 
events, the mixed events, and fitting the charged particle multiplicity distribution to the NBD, 
respectively. It is turned out that the correlation strength increases with increasing centrality bin 
size monotonously. If the discrepancy between real (total) and statistical correlation strengths is 
identified as dynamical one, the dynamical correlation may just be a few percent of the total (real) 
correlation. 
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I. INTRODUCTION 

The study of fluctuations and correlations has been 
suggested as a useful means for revealing the mechanism 
of particle production and Quark-Gluon-Plasma (QGP) 
formation in Relativistic Heavy Ion Collisions [liQ- Cor- 
relations and fluctuations of the thermodynamic quan- 
tities and/or the produced particle distributions may be 
significantly altered when the system undergoes phase 
transition from hadronic matter to quark-gluon matter 
because the degrees of freedom in two matters is very 
different. 

The experimental study of fluctuations and correla- 
tions becomes a hot topic in relativistic heavy ion colli- 
sions with the availability of high multiplicity event-by- 
event measurements at the CERN-SPS and BNL-RHIC 
experiments. An abundant experimental data have been 
reported [1, 0, [E] where a lot of new physics arise and are 
urgent to be studied. A lot of theoretical investigations 
have been reported as well [1, 0, S, i, [10, HH, 111 • 

Recently STAR collaboration have measured the 
charged particle forward-backward multiplicity correla- 
tion strength h in Au-|-Au collisions at ^snn=200 GeV 
[ll,[l3]- The outstanding features of STAR data are: 

• In most central collisions, the correlation strength 
b is approximately flat across a wide range in Ajy 
which is the distance between the centers of for- 
ward and backward (pseudo)rapidity bins. 

• This trend disappears slowly with decreasing cen- 
trality and approaches a exponential function of 
At; at the peripheral collisions. 

That has stimulated a lot of theoretical interests [1, 0, 
Hill. 
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In Ref. [lO[, a statistical model was proposed to cal- 
culate the charged particle forward-backward multiplic- 
ity correlation strength b in 0-10% most central Au+Au 
collisions at -ysNN=200 GeV. One outstanding feature 
of STAR data, the 6 as a function of A77 is approxi- 
mately flat, was well reproduced. The calculated value 
oib ^ 0.44 was compared with STAR data of 6 « 0.60 

EH. 

However, in this statistical model [l3| the Negative Bi- 
nomial Distribution (NBD) is assumed for the charged 
multiplicity distribution and the NBD parameters of /i 
and k (see later) are extracted from fit in with PHENIX 
charged particle multiplicity distribution [ll]. It is 
turned out in Ref. |ll| that the experimental 77 and 
Pt acceptances have large infiuences on the correlation 
strength b. The STAR experimental acceptances are 
quite different from PHENIX, thus the inconsistency, us- 
ing PHENIX multiplicit y d ata to explain STAR corre- 
lation data, involved in have to be studies further. 
Meanwhile, what is the discrepancy between b « 0.60 
(STAR datum) and 6 « 0.44 (NBD) also needs to be 
answered. 

In this paper we use a parton and hadron cascade 
model PACIAE fisj , to investigate the centrality bin size 
dependence of charged particle multiplicity correlation 
in 5, 0-5, and 0-10% most central Au-|-Au collisions at 
\/*NN=200 GeV. Following Ref. [ll| we generate the 
real events (6000) by the PACIAE model, construct the 
mixed events according to real events one by one, and 
extract the NBD parameters and k) from fitting the 
real events charged particle multiplicity distribution to 
the NBD. Then the charged particle forward-backward 
multiplicity correlation strength b is calculated for the 
real events (real correlation strength), the mixed events 
(statistical correlation strength), and the NBD (NBD 
correlation strength), respectively. They are all nearly 
flat across a widw range in Ar;. Their magnitude in 0- 
10% most central Au-f Au collisions are about 0.63, 0.59, 
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and 0.52, respectively. So the corresponding STAR data 
is well reproduced. It is turned out that the real (total), 
statistical, and NBD correlation strengths increase with 
increasing centrality bin size monotonously. If the dis- 
crepancy between real (total) and statistical correlation 
strengths is identified as dynamical one, then the dy- 
namical correlation strength may just be a few percent 
of the total (real) correlation strength. 

II. PACIAE MODEL 



collision assumed to be 40 mb. The corresponding col- 
lision time of this collision pair is then calculated. So 
the particle list and the NN collision (time) list can be 
constructed. A NN collision pair with smallest collision 
time is selected from the NN collision (time ) list and 
performed by the method in former paragraph. After 
upgrading the particle list and collision (time) list we se- 
lect and perform a new NN collision pair again. Repeat 
these processes until the collision (time) list is empty 
we obtain a initial partonic state for a nucleus-nucleus 
collision. 



The parton and hadron cascade model, PACIAE [16|, 
is based on PYTHIA [17] which is a model for hadron- 
hadron (hh) collisions. The PACIAE model is com- 
posed of four stages: parton initialization, parton evolu- 
tion (rescattering), hadronization, and hadron evolution 
(rescattering). 



A. Parton initialization 

In the PACIAE model a nuclcon-nuclcon (NN) col- 
lision is described with PYTHIA model, where a NN 
(hh) collision is decomposed into the parton-parton col- 
lisions. The hard parton-parton collision is described 
by the lowest-leading-order (LO) pQCD parton-parton 
cross section pjj with modification of parton distribu- 
tion function in the nucleon. And the soft parton-parton 
interaction is considered empirically. The semihard, be- 
tween hard and soft, QCD 2^2 processes are also 
involved in PYTHIA (PACIAE) model. Because of the 
initial- and final-state QCD radiation added to the above 
processes, the PYTHIA (PACIAE) model generates a 
multijet event for a NN [hh) collision. That is followed, 
in the PYTHIA model, by the string-based fragmenta- 
tion scheme (Lund model and/or Independent Fragmen- 
tation model) , thus a hadronic state is reached for a NN 
[hh] collision. However, in the PACIAE model above 
string fragmentation is switched off temporarily, so the 
result is a multijet event (composed of quark pairs, di- 
quark pairs and gluons) instead of a hadronic state. If 
the diquarks (anti-diquarks) are split forcibly into quarks 
(anti-quarks) randomly, the consequence of a NN (hh) 
collision is its initial partonic state composed of quarks, 
anti-quarks, and gluons. 

A nucleus-nucleus collision, in the PACIAE model, is 
decomposed into the nucleon-nuclcon collisions based on 
the collision geometry. A nucleon in the colliding nucleus 
is randomly distributed in the spatial coordinate space 
according to the Woods-Saxon distribution (r) and the 
47r uniform distribution [9 and 0). The beam momen- 
tum is given to Pz and px = Py — Q '\s assumed for each 
nucleon in the colliding nucleus. A closest approaching 
distance of two assumed straight line trajectories is cal- 
culated for each NN pair. If this distance is less than 
or equal to \/ CTtot/'"', then it is considered as a collision 
pair. Here atot refers to the total cross section of NN 



B. Parton evolution (rescattering) 

The next step, in the PACIAE model, is parton evolu- 
tion (partonic rescatterin g). Here the 2^2 LO-pQCD 
differential cross sections [18] are employed. The differ- 
ential cross section for a subprocess ij — > kl reads 



ij — 'kl 



dt 



= K- 



ij^kl 



(1) 



where the K factor is introduced for higher order correc- 
tions and the non-perturbative QCD correction as usual. 
Take the process qiq2 (71(72 as an example, one has 



E 



which can be regularized as 
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(2) 



(3) 



by introducing the parton colour screen mass, m=0.63 
GeV. In above equation s, t, and u are the Mandelstam 
variables and a.s= 0.47 stands for the running coupling 
constant. The total cross section of the parton collision 
i + i is 
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k,l 



dt- ^ 



dt 



(4) 



With these total and differential cross sections the par- 
ton evolution (rescattering) can be simulated by the 
Monte Carlo method until the parton-parton collision 
is ceased (partonic freeze-out). 



C. Hadronization 

In the PACIAE model the partons can be hadronized 
with the string-based fragmentation scheme or by the 
coalescence (recombination) models [H, [13, [HI, [2^, ■ 
The Lund string fragmentation regime, involved in the 
PYTHIA model, is adopted for hadronization in this pa- 
per, see p/7|] for the details. 
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Meanwhile, we have proposed a simulant coalescence 
(recombination) model which can be briefly explained as 
follows: 

1. The Field- Feynman parton generation mechanism 
[23 | is first applied to deexcite the energetic parton 
and thus to increase the parton multiplicity. This 
deexcitation of an energetic parton plays a similar 
role as string multiple fragmentation in the Lund 
model ,123]. 

2. The gluons are forcibly split into qq pair randomly. 

3. In the program there is a hadron table composed of 
mesons and baryons. The pseudoscalar and vector 
mesons made of u, d, s, and c quarks, as well as 
S+, 5°, B*°, and T are considered. The SU(4) 
multiplets of baryons made of u, d, s, and c quarks 
(except those with double c quarks) as well as A° 
are considered. 

4. Two partons can coalesce into a meson and three 
partons into a baryon (antibaryon) according to 
the flavor, momentum, and spatial coordinates of 
partons and the valence quark structure of hadron. 

5. When the coalescing partons can form either a 
pseudoscalar meson or a vector meson (e. g. ud 
can form either a tt"*" or a p"*") a judgment of less 
discrepancy between the invariant mass of coalesc- 
ing partons and the mass of coalesced hadron is in- 
voked to select one from two mesons above. In the 
case of baryon, e. g. both p and A+ are composed 
of uud, the same judgment is invoked to select one 
baryon from both of ^ and | baryons. 

6. Four momentum conservation is required. 

7. There is a phase space condition 



equal to y cr'^ot/^ l26|, then it is considered as a coUision 

pair. Here cr^o'j: refers to the total cross section of hh 
collision. The corresponding collision time of this colli- 
sion pair is then calculated. So the hh collision (time) 
list can be constructed. A hh collision pair with small- 
est collision time is selected from the collision (time) list 
and performed by the usual two-body collision method 
[2^ . After upgrading the particle list and collision (time) 
list we select and perform a new hh collision pair again. 
Repeat these processes until the collision (time) list is 
empty (hadronic freeze-out). 

A isospin averaged parametrization formula is used for 
the hh cross section [2^. However, we also provide 
a option of constant total, elastic, and inelastic cross 
sections ^tof = 40 mb, <7 = 25 mb, a^;^ = 

35 mb, crfj[ — 10 mb, and the assumed ratio of inelastic 
to total cross section equals 0.85. We also assume 



jpp 



= aAN = (tAA. 



(6) 



The cross section of ttN and kN, for instance, is as- 
sumed to be equal to the cross section of ttN and kN, 
respectively. 

The momentum of scattered particles in a hh elastic 
collision is simulated according to that the hh differential 
cross section, da^l[/dt, is assumed to be an exponential 
function of t which is squared momentum transfer [2^ . 
As it is impossible to include all inelastic channels, we 
consider only a part of them (« 600) which have notice- 
able effects on the hadronic final state, and the rest is 
attributed to the elastic scattering. Take incident chan- 
nel ttN as an example, if there are possible final channels 
of ttN ttA, TriV pN, and ttN fcA, their relative 
probabilities are then used to select one among above 
three channels. The momentum of scattered particles 
in a hh inelastic collision is simulated according to the 
usual two-body kinematics [2^ [2^ . 



(5) 



where h^/d is the volume occupied by a single 
hadron in the phase space, d=4 refers to the spin 
and parity degeneracies, Ar and Ap stand for the 
spatial and momentum distances between coalesc- 
ing partons, respectively. 



III. CALCULATION AND RESULT 

Following [s^l the charged particle forward-backward 
multiplicity correlation strength b is defined as 



(n/rib) - {nf){nb) _ cov(n/,n6) 



var(n/) 



(7) 



D. Hadron evolution (rescattering) 

We obtain a configuration of hadrons in spatial and 
momentum coordinate spaces for a nucleus-nucleus col- 
lision after the hadronization. If one only considers 
the rescattering among tt, k,p, n, p{uj), A, A, E, S, fi, J/^ 
and their antiparticles, the particle list is then con- 
structed by the above hadrons. A closest approaching 
distance of two assumed straight line trajectories is cal- 
culated for each hh pair. If this distance is less than or 



where Uf and Ub are, respectively, the number of charged 
particles in forward and backward pseudorapidity bins 
defined relatively and symmetrically to a given pseudora- 
pidity rj. {11 f) refers to the mean value of ny for instance. 
cov{nf ,nb) and var(n/) are the forward-backward multi- 
plicity covariance and forward multiplicity variance, re- 
spectively. 

In the calculations the default values given in the 
PYTHIA model are adopted for all model parameters 
except the parameters K and bg (in the Lund string 
fragmentation function). The ^=3 is assumed and the 
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TABLE I: Total charged particle multiplicity in three 77 
fiducial ranges in 0-6% most central Au+Au collisions at 
ysi^i7=200 GeV. 



NM < 4.7) < 5.4) iVch (total) 



PHOBOS" 
PACIAE 



4810 ± 240 
4819 



4960 ± 250 
4983 



5060 ± 250 
5100 



The experimental data are taken from [31] 



6s =6 is fixed by fitting the charged particle multiplicity 
to the corresponding PHOBOS data in 0-6% most cen- 
tral Au+Au coUisions at ^snn=200 GeV |31j as shown 
in Tab.lll Therefore in the generation of real events there 
is no free parameters. The mapping relation [3^ between 
the centrality definition in theory and experiment 



Ra + RB; 



(8) 



is employed. In the above equation bi (in fm) refers to 
the theoretical impact parameter and g stands for the 
percentage of geometrical (total) cross section used in 
experiment to define the centrality. Ra = 1.12 A^^^ + 
0.45 fm, for instance, is the radius of nucleus A. Thus 
the 0-10, 0-6, 0-5, and 5% most central collisions, for 
instance, are mapped to < 5^ < 4.46, < 6^ < 3.53, 
< 6i < 3.20, and bi = 3.20 fm, respectively. 

In this paper we propose a mixed event method where 
the mixed events are generated according to real events 
one by one. We first assume the charged particle multi- 
plicity n in a mixed event is the same as one correspond- 
ing real event. However, n particles of this mixed event 
are sampled randomly from the particle reservoir com- 
posed of all particles in all real events. Therefore, there 
is no dynamical relevance among the particles in a mixed 
event. So the correlation calculated by mixed events is 
reasonably to be identified as the statistical correlation 

m- 

It is known that the statistical correlation can also 
be studied by the NBD method, because the charged 
particle multiplicity distribution in high energy heavy- 
ion collisions is close to NBD [l5|. For an integer n the 
NBD reads 



P{n;fj.,k) 




(1 + /VA:)"+'= ' 



(9) 



where /x = (n) is a parameter, P{n; /i, k) is normalized 
in < n < 00, and k is another parameter responsible 
for the shape of the distribution. As proved in [llj] the 
correlation strength can be expressed as 



b = 



{nf)+k' 



(10) 



where the parameter k is fixed by fitting the charged 
particle multiplicity to the NBD usually. 

We compare the theoretical charged particle pseudo- 
rapidity distribution (open circles) in 0-6% most central 
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FIG. 1: Charged particle pseudorapidity distribution in 
Au+Au collision at y^SNN=200 GeV: (a) 0-6% most central 
collision and (b) 0-10, 0-5, and 5% most central collision. The 
experimental data are taken from j3ll ]. 



Au+Au collisions at ysNN=200 GeV with the corre- 
sponding PHOBOS data (solid squares) [3l| in Fig. 1 (a). 
One sees here that the PHOBOS data are well repro- 
duced. In Fig. 1 (b), we compare the charged particle 
pseudorapidity distributions in 0-5% (open circles) and 
5% (open triangles) most central Au+Au collisions with 
the 0-10% one (open squares). We see in Fig. 1 (b) that 
the pseudorapidity distribution in 5% most central col- 
lision is quite close to the 0-10% one, because the 5% 
centrality is nearly equal to the average centrality of 0- 
10% centrality bin. 

In Fig. 2 we compare the calculated real (total) cor- 
relation strength b (open squares) as a function of A77 
in 0-10% most central Au+Au collisions at ^snn=200 
GeV with the corresponding STAR data (solid squares) 
[13. The STAR data feature of correlation strength b 
is approximately flat across a wide range in Ar] are well 
reproduced. For comparison we also give the real (total) 
correlation strength in 0-5 and 5% most central collisions 
by open circles and triangles, respectively. One sees here 
that the real (total) correlation strength decreases with 
decreasing centrality bin size monotonously, because the 
charged particle multiplicity fluctuation decreases from 
0-10 to 0-5 and to 5% monotonously, as one will see in 
Fig. 3. This first result of the correlation strength in- 
creases with increasing centrality bin size monotonously 
given in the transport model remains to be proved ex- 
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FIG. 2: Charged particle forward-backward multiplicity cor- 
relation strength b in 0-10, 0-5, and 5% most central Au-|-Au 
collisions at a/snn=200 GeV. The experimental data are 
taken from [1^ . 
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FIG. 4: The calculated charged particle total (real), statisti- 
cal, and NBD correlation strengths in 0-10, 0-5 and, 5% most 
central Au+Au collisions at a/^nn=200 GeV. 
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FIG. 3: Charged particle multiplicity distributions in 0-10 
(open squares), 0-5 (open circles), and 5% (open triangles) 
most central Au+Au collisions at ysNN=200 GeV. The dot- 
ted, dashed, and solid lines are the corresponding NBD fits, 
respectively. 

The calculated charged particle multiplicity distribu- 
tions in 0-10, 0-5, and 5% most central Au+Au colli- 
sions at Y^SNN=200 GeV are given in Fig. 3, respec- 
tively, by the open squares, circles and triangles. The 
corresponding NBD fits are shown by dotted, dashed, 
and solid lines, respectively. One sees in Fig. 3 that the 
charged particle multiplicity fluctuation is increased and 
the NBD fit is worsened with increasing centrality bin 
size monotonously. 



In Fig. 4 we compare the calculated charged parti- 
cle real (solid symbols), statistical (open symbols), and 
NBD (lines) correlation strengths as a function of A77 
in 0-10, 0-5, and 5% most central Au+Au collisions at 
V^NN=200 GeV. The solid squares, open squares, and 
dotted line are for 0-10% most central collisions, solid 
circles, open circles, and dashed line for 0-5%, and solid 
triangles, open triangles, and solid line for 5%, respec- 
tively. We see in Fig. 4 that the behavior of correla- 
tion strength increases with increasing centrality bin size 
monotonously is not only existed in the real correlation 
strength but also in the statistical and NBD ones. 

If the discrepancy between real (total) and statisti- 
cal correlation strengths is identified as the dynamical 
correlation strength, one then sees in Fig. 4 that the dy- 
namical correlation strength may just be a few percent of 
the total (real) correlation strength. The dynamical cor- 
relation strength in 0-10% most central collision is close 
to the one in 5% most central collision globally speaking. 
That is because the later centrality is nearly the average 
of the former one. The dynamical correlation strength in 
0-10% most central collisions is globally less than 0-5% 
most central collision. That is because the interactions 
(represented by the collision number for instance) in the 
former collisions is weaker than the later one. We also 
see in Fig. 4 that the statistical correlation strength is 
nearly the same as the NBD one in the 5% most cen- 
tral collision, that is consistent with the results in p + p 
collisions at the same energy [11] . However the discrep- 
ancy between statistical and NBD correlation strengths 
seems to be increased with increasing centrality bin size 
monotonously. That is mainly because the NBD fitting 
to the charged particle multiplicity distribution becomes 
worse with increasing centrality bin size monotonously. 
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IV. CONCLUSION 

In summary, we have used a parton and hadron 
cascade model, PACIAE, to study the centraUty bin 
size dependence of charged particle forward-backward 
multiplicity correlation strength in 5, 0-5, and 0-10% 
most central Au-f Au collisions at -^snn^SOO GeV. The 
real (total), statistical, and NBD correlation strengths 
are calculated by real events, mixed events, and NBD 
method, respectively. The corresponding STAR data 
feature of the correlation strength b is approximately flat 
across a wide range in Arj in most central Au+Au col- 
lisions is well reproduced. It is turned out that the cor- 
relation strength increases with increasing centrality bin 
size monotonously. This first result, given in the trans- 



port model, remains to be proved experimentally. If the 
discrepancy between real (total) and statistical correla- 
tion strengths is identified as dynamical one [ll|, then 
the dynamical correlation may be just a few percent of 
the total (real) correlation. As a next step, we will in- 
vestigate the relation between correlation strength b and 
the centrality bin size in the mid-central and peripheral 
collisions, and the STAR data feature of b approaches an 
exponential function of Arj at the peripheral collisions. 
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